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Abstract 

Nonlinear z-independent solutions to a generalized Grad-Shafranov 
equation (GSE) with up to quartic flux terms in the free functions 
and incompressible plasma flow non parallel to the magnetic field are 
constructed quasi-analytically. Through an ansatz the GSE is trans- 
formed to a set of three ordinary differential equations and a constraint 
for three functions of the coordinate x, in cartesian coordinates (x, y), 
which then are solved numerically. Equilibrium configurations for cer- 
tain values of the integration constants are displayed. Examination of 
their characteristics in connection with the impact of nonlinearity and 
sheared flow indicates that these equilibria are consistent with the L-H 
transition phenomenology. For flows parallel to the magnetic field one 
equilibrium corresponding to the H-state is potentially stable in the 
sense that a sufficient condition for linear stability is satisfied in an 
appreciable part of the plasma while another solution corresponding 
to the L-state does not satisfy the condition. The results indicate that 
the sheared flow in conjunction with the equilibrium nonlinearity play 
a stabilizing role. 
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1 Introduction 



For axisymmetric toroidal plasma equilibria the force-balance equation and 
Maxwell's equations reduce to the Grad-Shafranov equation (GSE) for the 
poloidal magnetic flux function ip [TJ, [2]. Analytical solutions to the GSE 
are obtained by specifying the plasma pressure and poloidal current functions 
of usually in such a way as to linearize the resulting partial differential 
equation, e.g. [3J-[T2]. Analytical solutions to the GSE are very useful for 
theoretical studies of plasma equilibrium, transport and stability as well as 
benchmarks for numerical codes [T3] - Also it has been established in a va- 
riety of magnetic configurations that sheared flows can reduce turbulence 
and produce transport barriers, which under certain conditions can extend 
to the whole plasma, e.g. [14]. In view of a fusion reactor the spontaneous 
formation of transport barriers, i.e., those driven by internal processes even 
in the absence of external sources, is of particular interest. For this reason 
among others stationary equilibria with plasma flow have been extensively 
studied on the basis of generalized GSEs, e.g. [H)]-[37]. In particular, al- 
though complex numerical codes are extensively used to attempt simulations 
of the L-H transition, certain equilibrium considerations in connection with 
this transition are helpful, e.g. [38]-|41j. 

The simplest known and widely used in various studies, analytical solu- 
tion to the GSE, is the Solovev equilibrium [3J. Extension of the original 
Solovev solution, to include the possibility of sheared flows appeared in |23j. 
In other extensions additional free parameters were introduced to construct 
configurations with fusion relevant plasma boundaries and desirable values of 
confinement figures of merit as the safety factor on magnetic axis [TTJ 021 US] • 
Most of the solutions are associated with pressure and current profiles, in- 
cluding up to quadratic terms in the flux function ip to linearize the resulting 
equation [3J-[T2]. Linear equilibria with flow were constructed in |15]-j35] and 
in Refs. cited therein. Also, the non linear translational symmetric equilibria 
of "cat eyes" and counter rotating vortices with flow parallel to the magnetic 
field were studied in [3SJ |37]. These nonlinear equilibria, however, are peri- 
odic in one direction (x) and therefore the plasma is not bounded along this 
direction. 

In most of the above cases the axisymmetric equilibria are obtained as 
separable solutions of GSE. A novel non-separable class of solutions was 
found in Ref. [30] describing up-down symmetric configurations with incom- 
pressible flows parallel to the magnetic field and it was extended recently to 
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include asymmetric configurations [21] and flows of arbitrary direction [35J. 
For non parallel flows the question of the stability is usually not considered 
and this is partly due to the difficulty of the subject and the absence of a 
concise criterion. Few sufficient conditions for linear stability are available 
only for parallel flows [42j-[44j. In previous studies we found that the sta- 
bility condition of Ref. [33] is not satisfied for the linear equilibria of [27] 
and [34J while it is satisfied within an appreciable part of the plasma for 
the nonlinear equilibria of |36j and [37] . This led us to the conjecture that 
the equilibrium nonlinearity may act synergetically with the sheared flow to 
stabilize the plasma. 

Aim of the present study is to construct certain two dimensional non- 
linear up-down symmetric equilibria with incompressible flow of arbitrary 
direction in z-independent geometry. They are more pertinent to a magneti- 
cally confined plasma than those of Refs. [36J and [37] because the plasma is 
bounded on the poloidal plane. Another reason for considering translational 
symmetry is the many free physical and geometrical parameters involved in 
connection with the flow amplitude, direction and shear, equilibrium non- 
linearity, symmetry and toroidicity. Thus, in the presence of nonlinearity 
one first could exclude toroidicity. The study is performed quasi-analytically 
through a non separable ansatz under which the GSE is transformed to a 
set of three ordinary differential equations and a constraint for three func- 
tions. The solutions give nested magnetic surfaces and their characteristics 
are studied by means of certain equilibrium quantities and confinement fig- 
ures of merit as the safety factor, electric field and E x B velocity shear. Also, 
for parallel flows the linear stability is considered by means of the aforemen- 
tioned sufficient condition [33] . The results are in qualitative agreement with 
phenomenological characteristics of an edge transport barrier, confirm rele- 
vant scenarios [13], [23J and support the above conjecture. 

The organization of the paper is as follows: In the first section we briefly 
review the general setting for the equations of incompressible flow with trans- 
lational symmetry together with the generalized GSE. In Section II the pro- 
posed ansatz and the resulting equations are presented which then are inte- 
grated numerically. In section III we consider the solutions for certain values 
of the various parameters and integration constants and discuss the most 
important equilibrium properties. In section IV the criterion for linear sta- 
bility is applied to the equilibria constructed for parallel flows. Section V 
summarizes the study and briefly proposes potential extensions. 
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2 Translational Symmetric Equilibria with flow 



The equilibrium of a cylindrical plasma with incompressible flow and 
arbitrary cross-sectional shape, satisfies [H], [23J, 

(1 - M p 2 )V 2 ^ - \{Ml)'\V^ + A ^ oPs + ^pj=0 (1) 

for the poloidal magnetic flux function if). Here, M p (if>), P s (ij)), p(if)) and 
B z (if>) are respectively the poloidal Alfven Mach function, pressure in the 
absence of flow, density and magnetic field parallel to the symmetry axis 
z, which are surface quantities. Because of the symmetry, the equilibrium 
quantities are z— independent and the axial velocity v z does not appear ex- 
plicitly in (TfJ). Derivation of ([T]) is based on the following two steps: First, 
express the divergence free fields in terms of scalar quantities as 

B = B z Vz + Vz xVip 
Pol = V 2 ipVz -Vzx VB Z 
pv = pv z Vz + Vz x VF 

and the electric field by E = — V$. Second, project the momentum equation, 
p (v • V) v = j x B — VP, and Ohm's law, E + v x B = 0, along the symmetry 
direction z, B and Vip. The projections yield four first integrals in the form 
of surface quantities (two out of which are F{ip) and $(^)), Eq. flTJ and the 
Bernoulli relation for the pressure 

p = p s (ip) - ^-Mim^\ 2 (2) 

Because of the flow P is not a surface quantity. Also the density becomes 
surface quantity because of incompressibility and Mp(if>) = (F (ip)) 2 / (pop) ■ 
Five of the surface quantities, chosen here to be P s , p, B z , M 2 and v z , 
remain arbitrary. 

Using the transformation 

= J\l - Mfrgftftdg, {Ml < 1) (3) 
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Eq. flTJ is transformed to 

Note that transformation ([3]) does not affact the magnetic surfaces, it just re- 
labels them. Eq. fll]) is identical in form with the static equilibrium equation. 
In the present study we assign the free function term in (J3J) as 

( o B l\ U " U " ^ 

I H0 P s + — I = C + CiM + C 2 y + C 3 y + C 4 — (5) 

where Co, ci, . . . , C4 are free parameters. 

3 Proposed Ansatz 

We use Eq. (jSJ) into Eq. (@|, employ the ansatz 

iy^ly' + ZM-PoM ... 

"= y^W) <6) 

and equate the nominator of the resulting equation to zero. From the 
y 6 — terms we obtain (a prime denotes derivative with respect to x) 

N'l + a + CaiVx + c 3 Nl + c 4 Nl = (7) 

From the y°— terms we obtain the constraint C s = 0, where 

C s = 2(Nx -f) + D [ Cl + c 2 f + c 3 f 2 + c 4 / 3 ] = (8) 

The y A and y 2 — terms, after rearrangement yield 



(9) 



and 



V2 (iW) V2 X" 6 + 

+c 3 D (N 1 - /) + 3c 4 D iVi(iVi - /) = (10) 
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Eq. ([7[) is solved using the tanh method |45| . a method of solving non linear 
differential equations, which also employed in [2S]. We have two solutions. 
The first is N\(x) = ao + axtanh(vx), where 

ci + c 2 a + c 3 aQ + c 4 ao = 

c 2 + c 3 (2a ) + c 4 (3ao) = 2v 2 

c 3 + c 4 (3a ) = 

c 4 (a 2 ) = -2v 2 (11) 

and the second is Ni(x) = a + ai/cosh(vx), where 

ci + c 2 a + c 3 aQ + c 4 al = 

c 2 + c 3 (2a ) + c 4 (3ag) = -v 2 

c 3 + c 4 (3a ) = 

c 4 (a 2 ) = 2v 2 (12) 



4 Solutions and equilibrium properties 



We have solved numerically Eqs. flS}, © and f flQl) . Using the first of the 
solutions for Ni, namely the tanh solution, we obtained the equilibrium of 
Fig. 1. We have used a = 1.1, ai = 2.5, v = 0.6 and in Eq. (JSJ c = 2.588, 
Ci = —0.638, c 2 = 0.302, c 3 = 0.38, c 4 = —0.115. The boundary flux surface 
corresponds to = 0.11 while on the magnetic axis u a = 0. The constraint 
was kept close to zero for the whole of the integration process and we got an 
average value of \C S \ equal to 0.10. Given the nonlinearity and complexity 
of the method this implies that the solution is indeed acceptable. Simple 
quadratic fitting gives 

/ = 1.272z 2 + 0.049x + 0.001 and D = 1.488a: + 3.3. 

Using the second of the solutions for Ni, namely the cosh solution, we 
obtained the equilibrium of Fig. 2. We have used ao = 1.0, a 4 = —1.6, 
v = 1.15 and in Eq. © c = 2.588, c x = 0.289, c 2 = 1.777, c 3 = -3.099, 
c 4 = 1.033. The boundary flux surface corresponds to u b = —0.05 while on 
the magnetic axis u a = 0. The constraint was kept close to zero for the whole 
of the integration process and the average value of |C S | was 0.01. Simple 
quadratic fitting gives / = -0.542x 2 + 0.009x and D = 0.994x 2 + 3.3. 
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Here instead of the velocity v z we have used the axial Mach function, 



Ml 



u 



v 2 / (B 2 /(fi p), and the approximation Mj 



M 2 
p 



relation to the tokamak scaling B„ « 0.1B Z and v„ ~ 0.1v 2 



In addition, to 



completely construct the equilibrium we have made the following choices 
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(13) 
(14) 

(15) 
(16) 



for the poloidal Mach function, axial Mach function, axial magnetic field and 
density, respectively, with B z0 = 2.24 T, p a = 4 x 10~ 7 Kgr/m 3 7 = 0.02, 
u a = 0, u b = 0.11 W m 2 , (with the subscripts a and b indicating the magnetic 
axis and boundary respectively), A = 0.5, m = 9n and M za = 1.1 M„„. with 



L pa 

various values of the parameters M pa and n. Here, Eqs. ffl3|) and fTl4|) can 
decribe Mach functions localized in the edge plasma region in connection 
with the L-H transition (in particular flows localized nearly in the one tenth 
of the exterior plasma will be considerded as it is shown in Fig. 3); Eq. ( fl~5|) 
represents a diamagnetic B z (u) (Fig. 4). Then, ([2]) and ()5]) imply a pressure 
peaked on axis (Fig. 5). 

Furthermore we have examined certain equilibrium characteristics by 
means of the safety factor, magnetic shear, axial current density, radial elec- 
tric field and E x B velocity shear, , and found the following results. 

1. The safety factor for both solutions shown in Figs. 6,7 is slightly af- 
fected by the flow. Also, the flow affects slightly the magnetic shear 
given by s(u) = 2(V/q)(dq/dV) as it can be seen in Fig. 8 for Equilib- 
rium 1. A similar plot holds for Equilibrium 2. 

2. The radial electric field for the two solutions has an extremum in the 
edge region which increases with flow (Figs. 9 and 10). The position of 
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the extremum, however, is nearly unaffected by the flow. These char- 
acteristics are indicative that the solutions may be relevant to the L-H 
transition as discussed in [23] where a similar behavior of the electric 
field was found (Fig. 3 therein). 

3. The E x B velocity shear which is believed to play a role in the tran- 
sitions to improved confinement regimes of magnetically confined plas- 
mas is given by 



d 


E x B" 




dr 


B 2 





(17) 



where r is the length variable normal to the magnetic surfaces. For 
Equilibrium 1 it is plotted in Fig. 11; a similar plot holds for Equilib- 
rium 2. coexb is increased by the flow in the edge region outer from the 
local minimum while it remains nearly unaffected in the central region. 
This is another indication supporting the relevance of the solutions to 
the L-H transition. 

4. The flow makes the axial ("toroidal") current density profile hollow as 
shown in Fig. 12 for Equilibrium 1. (A similar j tor profile is found for 
Equilibrium 2.) The larger the flow is the stronger the hollowness. Hol- 
low jtor profiles are usually related to the formation of internal transport 
barriers in tokamaks. However, despite of this characteristic and the 
fact that ojexb becomes maximum on the magnetic axis (Fig. 11) it 
is unlikely that the present equilibria are related to internal transport 
barriers because the safety factor is monotonically increasing from the 
magnetic axis to the plasma edge (Figs. 6, 7). According to observa- 
tions in tokamaks, e.g jl6] for JET and (17] for DIII-D, it is the reversed 
magnetic shear which plays a role in triggering the ITBs development. 
Also, as can be seen in Figs. 6 and 7 the flow makes the central q- values 
lower. 

5 Stability consideration 

We now consider the important issue of the stability of the solutions 
constructed in Section IV with respect to small linear MHD perturbations 
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by applying the sufficient condition of Ref. jS]. This condition states that a 
general steady state of a plasma of constant density and incompressible flow 
parallel to B is linearly stable to small three-dimensional perturbations if the 
flow is sub-Alfvenic (M 2 < 1) and A > 0, where A is given below by ffl8|) . 
Consequently, using henceforth dimensionless quantities we set p = 1. Also, 
for parallel flows (v = MB) it holds M p = M z = M. In fact if the density 
is uniform at equilibrium it remains so at the perturbed state because of 
incompressibility |48] . In the w-space for axisymmetric equilibria A assumes 
the form 

A = -f 



+ 

+ g- 

with 



j x V«) • (B ■ V)Vu+ 

| Vu- V(5 2 /2)+ 



'M 2 p \ \\/u\ 2 




(1 - M 2 ) 3 / 2 

ill 

1 - M 2 )V2 J 

P' s (u)-(M 2 )'B 2 /2 



(18) 
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1-M 2 



Symbolic computation of A over a wide rage of parametric values led to the 
following results: 

1. Equilibrium 1 is not satisfied, since A < everywhere, while Equilib- 
rium 2 is satisfied in an appreciable part of the plasma region. However, 
it is noted that since the stability condition is necessary, A < does not 
imply that an equilibrium is unstable. An example of the sign of A for 
Equilibrium 2 is given in the three-dimensional plot of Fig. 13. Also, 
profiles of A in the middle-plane y = for a static and a stationary 
equilibrium are shown in Fig. 14. 

2. Increase of M pa makes A more positive in the edge region as can be 
seen in the example of Fig. 14. A similar impact on A has the flow 
shear parameter n (Eq. (Fl3|) ) as can be seen in Fig. 15 showing the 
profile of 5 A = A(y = 0, n = 2) - A(y = 0, n = 1). 

3. The equilibrium nonlinearity in connection with the parameters C3 and 
C4 has a stabilizing effect in the edge region as shown in the example 
of Fig. 16 plotting the profile of the difference 5 A between a nonlinear 
and a linear Equilibrium 2. 
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According to the above results and the believe that the sheared flow is de- 
veloped during the L-H transition we conjecture that a static Equilibrium 1 
could correspond to the L state and a stationary Equilibrium 2 with ExB/5 2 
velocity shear to the H state. In a quasistatic evolution approximation the 
plasma could then evolve through successive states with increased sheared 
flow (increasing values of the parameters M pa and M za and most importantly 
increasing values of the shearing parameters m and n). 




-0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 

x 



Figure 1: Equilibrium 1. The bounding flux surface corresponds to Uf, = 0.11, 
with u a = 0, for the magnetic axis. For this equilibrium the average value of 
\C S \ is 0.10. 



6 Summary 

Two classes of solutions of nonlinear two dimensional magnetohydrodynamic 
equilibria for bounded magnetically confined plasmas with sheared incom- 



10 




-0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 



x 

Figure 2: Equilibrium 2. The bounding flux surface corresponds to Ub = 
—0.05, with u a = 0, for the magnetic axis. For this equilibrium the average 
value of \C S \ is 0.10. 

pressible non parallel flows have been constructed in cylindrical (z-independent) 
geometry. The equilibria hold for four arbitrary surface functions which were 
chosen to be the plasma density, axial Mach function, poloidal Mach-function 
and static pressure. 

After assigning the free functions, a systematic examination of equilib- 
rium quantities and confinement figures of merit, as the safety factor, electric 
field and E x B velocity shear for a variety of parametric values, implies that 
the equilibrium characteristics are qualitatively consistent with experimental 
evidence of the L-H transition. In addition, application of a sufficient condi- 
tion for linear stability and parallel flow indicates that one statioanry equi- 
librium being potential stable may describe the H-state and another static 
equilibrium not satisfying the stability condition the L-state. In addition 
the equilibrium non-linearity in conjunction with the flow and the flow shear 
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0.1 




x 



Figure 3: L-H transition-like Mach function in connection with Eq. ( {TBI with 
n = 1 and a maximum localized at a distance from the boundary as large as 
the on tenght of the minor radius. 

may play a stabilizing role. Although understanding the physics of the L-H 
transition remains incomplete the results of the present study may shed some 
light towards that goal. 

Finally it would be interesting trying to generalize these classes of solu- 
tions to up-down asymmetric configuration with a lower x- point in connec- 
tion with the ITER project. Also the study could be extended to toroidal 
geometry in order to examine the impact of toroidicity. 
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Figure 13: Stability function for the solution of Equilibrium 2. For most part 
of the equilibrium is positive and assumes negative values only in the core of 
the equilibrium. 
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Figure 14: Plot of the stability function A, at y=0, for the second equilib- 
rium for nonzero values of the nonlinearity parameters c3=-3.099, c4=1.033. 
Increasing the flow parameter M pa appears to improve stability for most part 
on the middle-plane except for the canter of the equilibrium 
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Figure 15: Plot of the deference S A = A(n = ...)— A(n = ...) for c3 = —3.099, 
c4 = 1.033 and Mpa = 0.1 clearly indicating that the stability is improved 
(SA > 0) at the external part of the middle-plane y = as the flow-shear 
parameter n increases. The dotted line represents the a;-axis. 
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Figure 16: Plot of the profile 6 A = A(c 3 = -3.099, c 4 = 1.033) - A(c 3 = c 4 
0) at y = indicating that the nonlinearity has a stabilizing effect (5 A > 
in the edge region 
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